Miscibility in a degenerate fermionic mixture induced by linear coupling 
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We consider a one-dimensional mean-field-hydrodynamic model of a two-component degener- 
ate Fermi gas in an external trap, each component representing a spin state of the same atom. 
We demonstrate that the interconversion between them (linear coupling), imposed by a resonant 
electromagnetic wave, transforms the immiscible binary gas into a miscible state, if the coupling 
constant, k, exceeds a critical value, K cr . The effect is predicted in a variational approximation, 
and confirmed by numerical solutions. Unlike the recently studied model of a binary BEC with 
the linear coupling, the components in the immiscible phase of the binary fermion mixture never 
fill two separated domains with a wall between them, but rather form anti-locked (7r-phase-shifted) 
density waves. Another difference from the bosonic mixture is spontaneous breaking of symmetry 
between the two components in terms of numbers of atoms in them, Ni and N2 . The latter effect is 
characterized by the parameter v = (Ni — N3)/ (Ni + N2) (only Ni + N2 is a conserved quantity), 
the onset of miscibility at n > k ci meaning a transition to v = 0. At k < n CI , v features damped 
oscillations as a function of k. We also briefly consider an asymmetric model, with a chemical- 
potential difference between the two components. The relation between the imbalance in the spin 
population, induced by the linear coupling, and the developing spatial structure resembles the known 
Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) states in the Fermi mixture. Dynamical states, when k 
is suddenly switched from zero to a value exceeding « cr , are considered too. In the latter case, the 
system features oscillatory relaxation to the mixed state. 

PACS numbers: 03.75.Ss 
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I. INTRODUCTION 

Experimental and theoretical studies of degenerate 
Bose [lj and Fermi [2| gases at ultra-low temperatures 
have become a vast research area. This field offers a 
unique testbed for the implementation of various ideas 
from classical and quantum nonlinear dynamics. In par- 
ticular, experiments are possible in binary mixtures of 
atoms in two different hyperfine (spin) states. The mix- 
tures were created by means of the sympathetic-cooling 
technique in the Bose-Einstein condensate (BEC) made 
of 87 Rb atoms In such a binary condensate, mu- 
tual interconversion, i.e., linear coupling between the two 
species, can be induced by means of a resonant (spin- 
flipping) electromagnetic wave, with a frequency in the 
GHz range 0]. Various effects were considered in this 
setting, including Josephson oscillations between the two 
states domain walls between them "breathing- 
together" oscillatory modes in the mixture 7], and non- 
topological vortices 

In a recent work , a new effect was predicted, namely, 
transition to miscibility induced by the linear coupling 
between the immiscible components in a binary conden- 
sate. It was demonstrated that, essentially, it occurs as 
a transition from a state with the two components sepa- 
rated by a domain wall to a completely mixed one. Quan- 
titatively, the transition is characterized by an order pa- 
rameter which determines the separation between centers 
of mass of the two components. With the increase of the 



linear-coupling strength (coupling constant), the order 
parameter decreases, and the full mixing sets in when it 
vanishes. Experimental observation of this effect is fea- 
sible due to the fact that the mixture of two hyperfine 
states in 87 Rb is immiscible, but close to the miscibility 
threshold 0, El ■ The proximity of the mixture to the 
threshold can be additionally controlled by means of the 
Feshbach resonance, which can be induced by an external 
magnetic field 11] or optical pump |T^ . 

Phase separation (demixing) due to the repulsion be- 
tween two components was also investigated in boson- 
fermion mixtures |13|. and, more recently, in a mixture 
of two fermion species [LH Ha ] . In the latter case, a situ- 
ation is possible when the species are two different spin 
states of the same fermion atom (as realized experimen- 
tally in 40 K and 6 Li [13), hence the linear intercon- 
version between them can be induced by the properly 
tuned electromagnetic wave. 

In a recent experimental work 0] , the population im- 
balance between two spin states ( "polarization" ) in a de- 
generate 6 Li gas (trapped in a cigar-shaped cavity) was 
controlled by means of this technique. It was observed 
that, due to strong attraction between the different states 
(induced by the Feshbach resonance in external magnetic 
field) they formed a fully mixed state in the central part 
of the trap, that was separated by two sharp boundaries 
from side domains, filled with the atoms representing a 
majority state (see Figs. 3 and 5 in Ref. [l8|). 

Thus, it is relevant to consider the transition to mixing 
in the binary fermionic gas under the action of the linear 
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coupling (as well as in the bosonic mixture, the proximity 
to the miscibility point may be additionally controlled 
by means of the Feshbach resonance, which affects the 
strength of the repulsion between the two species [T^). 
This is the subject of the present work. 

Rigorous description of the dynamics of the degener- 
ate gas of fermions is based on the system of quantum- 
mechanical equations of motion for individual fermions 
(see, e.g., Ref. j2(jj), which is difficult to handle if the 
gas contains many atoms. Simplified approaches, of the 
mean-field-hydrodynamic (MFHD) type, were also elab- 
orated for the description of quantum fermionic gases, 
in static [3 and dynamical [H IH H]| settings. In 
the latter case, the approximation amounts to a single 
evolution equation for the fermion wave function, which 
formally resembles the Gross-Pitaevskii equation (GPE) 
for the boson gas , but features a different nonlinearity. 
This MFHD equation was successfully used to describe 
collapse |22j and solitons |23] in degenerate fermionic 
gases, and a system of nonlinearly coupled equations was 
employed to predict the phase separation in the binary 
gas with repulsion between the species [l5| . 

In this work, we aim to investigate the onset of mix- 
ing in the immiscible binary fermion gas under the ac- 
tion of the extra linear coupling, in the (effectively) one- 
dimensional situation. In Sec. II, we formulate the 
model, and then elaborate an analytical approach to the 
problem, based on a variational approximation (VA). The 
latter provides a description of the binary mixture in 
terms of the number-of-atoms ratio in the two compo- 
nents, N1/N2. In the presence of the linear mixing, only 
the total number is conserved, N\ + N2, while the ratio 
v = (Ni — N2) I (Ni + N 2 ), if different from zero, char- 
acterizes spontaneous breaking of symmetry between the 
components. In fact, v ^ is a specific, and most impor- 
tant, manifestation of the immiscibility in the fermionic 
system with the linear coupling (in the usual setting, 
without the linear interconversion, N\ and N2 are con- 
served separately, which makes the situation altogether 
different). The VA predicts the onset of mixing at a 
critical value of the coupling constant, k = K cr , with 
v = at k > K cr . The critical value is predicted as a 
a function of parameters measuring the total number of 
atoms and strength of the nonlinear repulsion between 
the species. Comparison with numerical results demon- 
strates that VA provides a qualitatively correct descrip- 
tion of the transition to miscibility. 

Basic numerical results are presented in Sec. III. They 
demonstrate a drastic difference of the phase separa- 
tion in the fermionic mixture from the immiscible bi- 
nary BEC: while in the boson gas the two species spon- 
taneously separate and form a domain wall, with each 
component occupying the domain on either side of the 
wall (both in the ordinary setting, without the linear 
coupling [24|, and when the coupling is switched on, 
up to the onset of the induced mixing [9|). A binary 
fermionic gas features the separation in an altogether dif- 
ferent form: each component forms a spatially modulated 



symmetric (even) density wave, which are anti-locked, 
with the phase shift of tt between the components. A 
similar picture was predicted in the fermionic binary gas 
mixture without the linear coupling |l5j| . Here, a major 
novelty is the above-mentioned overall spontaneous sym- 
metry breaking between the components, manifested by 
N1/N2 7^ 1 (in the immiscible bosonic condensate, the 
linear coupling does not break the equality between Ni 
and N2 H ). As k approaches K cr (the mixed state), we 
observe N1/N2 — > 1, while the modulation depth in the 
density waves decreases. 

In Sec. IV, we briefly consider additional issues. One of 
them is the influence of asymmetry between the species, 
in the form of a chemical-potential difference between 
them, that can be induced by coupling of the atomic 
spin to external DC magnetic field. Also considered is 
a dynamical relaxation of the binary gas into a mixed 
state, when n is suddenly switched from zero to a value 
exceeding K cr . The results are summarized in Sec. V. 

To conclude the Introduction, it is relevant to men- 
tion that the relation between the imbalance in the spin 
population, induced by the linear coupling, and the spa- 
tial density-modulation waves developed by the two in- 
teracting species is similar to the known LOFF (Larkin- 
Ovchinnikov-Fulde-Ferrell) phase in a mixture of two 
fermion species with strongly differing values of the Fermi 
radius. As shown in original works |25j (see also review 
|26|). in that situation the pairing between the fermions 
belonging to the different species gives rise to Cooper 
pairs with a nonzero momentum, which translates into 
spatial structures in the form of crystalline patterns of 
the gap parameter. The LO FF p hase has various impli- 
cations for quark matter |26l |27| , and a possibility of its 
realization in the mixture of two different spin states in 
an ultracold fermion gas, with the interaction between 
them controlled via the Feshbach resonance, has been 
analyzed recently j2^| . 

II. THE MODEL AND VARIATIONAL 
APPROXIMATION 

A. The mean-field hydrodynamic (MFHD) 
equations 

Our starting point is the system of nonlinearly coupled 
one-dimensional MFHD equations for the fermion wave 
functions (f> and ip, as worked out in Ref. 15j, to which 
we add, following Ref. [9j and references therein, the 
linear-coupling terms with strength k and, in the general 
case, a chemical-potential difference, A/u (in the larger 
part of the paper, we consider the symmetric system, 
with Afj, = 0). In a scaled form, the equations are: 

%<k = -$ xx +x 2 (j)+(g l \(j)\ i l z + g 2 \ii\ 2 ^-^,{l) 

iipt = -i/j xx + x 2 tp+ (3iM 4/3 +S2H 2 ) tp-ncj) 
+ Afi-ip. (2) 
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Here, we imply the normalization 



N 



-Hoc. 



[<S> 2 (x) + ^ 2 (x)] dx = N 1 + N 2 



(3) 



Hence the total number of atoms is hidden in the posi- 
tive self-repulsion parameter |l5j| . gi, coefficient g 2 > 
accounts for the nonlinear repulsion between the the 
species, and the strength of the trapping potential (the 
coefficient in front of the parabolic potential, x 2 ) is scaled 
to be 1. The use of the same diagonal nonlinearity co- 
efficient gi in both equations implies the same number 
of atoms in them, and, unless we introduce k / 0, the 
number of atoms will remain equal in both channels. A 
characteristic feature of the MFHD approximation for 
the fermion gas is the power, 4/3, of the nonlinear self- 
repulsive terms (recall the power is 2 in the ordinary 
GPE). 

In the particular case of gi — 0, Eqs. and @ 
admit a different interpretation, being tantamount to a 
system of GPEs for the binary BEC [3 , in which the self- 
interaction coefficient (intra-species scattering length) is 
set equal to zero by means of a Feshbach-resonance ad- 
justment. This particular case establishes a link of the 
considered model to those for boson mixtures. 

Coefficient g 2 in Eqs. and ©, which accounts for 
the inter-species repulsion, can be controlled by means of 
the Feshbach-resonance technique in the fermion binary 
gas [l^. On the other hand, coefficient <?i accounts for 
the Fermi pressure and cannot be altered by means of 
this technique. In the analysis presented below, we vary 
both coefficients 171 and g 2 , which corresponds to the nor- 
malization adopted in Eqs. Q an d @ (in particular, 
varying gi actually implies the change in the number of 
atoms in the gas, as aid above). 

We look for stationary solutions, {ip, </>} — 
{^(x), &(x)} e~ HLt , with chemical potential n and 
real functions and $ obeying the equations 



M $ = -$" + a; 2 $+ (pi* 473 +.g 2 * 2 j k*, (4) 
= + x 2 ^ + (gi^ 4/3 + ,g 2 * 2 ) * - k$ 



(5) 



(the prime stands for d/dx). Together with the normal- 
ization condition, f_ °° [$ 2 (x) + \l/ 2 (x)] dx — 2, as per 
Eq. |J5J, Eqs. @ and (JSJ can be derived as the varia- 
tional equations, 5L/S& = 8L/5^ = dL/dfi = 0, from 
the Lagrangian, 



L 



x 2 ($ 2 + * 2 ) 



M ($ 2 + * 2 ) + (V) 2 + (*') 2 ) 



I 91 



$10/3 + ^10/3\ +ff2 $2^ 



2k$* + -A/i • W 



dx. 



(6) 



As said above, spontaneous symmetry breaking be- 
tween the two species in the system with the linear cou- 



pling, where norms N% and N 2 are not conserved sepa- 
rately, can be characterized by the ratio, 



iVi- A^ 2 
Ni +N 2 



[<J> 2 (z)-* 2 (z)]. (7) 



In Eq. Q, the normalization of Ni + N 2 is taken into 
account as per Eq. ©. The objective of the analysis 
is to find v as a function of k (for given g\ and g 2 ), and 
thus identify the above-mentioned critical value, K cr , past 
which (at k > K cr ) only the symmetric solution exists, 
with v = [the symmetric solution exists at k < n CI 
too, but it is unstable (it does not represent the system's 
ground state) in that case, cf. a similar situation in the 
binary BEC with the linear coupling 



B. Variational approximation 

The only tractable ansatz with which the VA can be 
applied to Lagrangian |JB} is based on simple Gaussians, 



$ = ae-* 2 / 2 , * = be-* 2 ' 2 , 



(8) 



with constants a and b. Normalization condition ||5J for 
this ansatz takes the form of 



(a 2 + b 2 



(9) 



and the symmetry-breaking parameter, defined by Eq. 
0, is 



V = (y/n/2) (a 2 



(10) 



Following the general method of the VA [23], we sub- 
stitute ansatz (jSJ in Lagrangian © and perform the in- 
tegration, which yields 



L 



(ab) — nab. 



ZVlbgi 
50 



2 v / 2 



(a 10 / 3 + 6 10 / 3 
(11) 



(normalization condition © can also be derived from 
this Lagrangian, as dL/d/j, = 0). In the framework of 
the VA, constants a and b are determined by equations 
dL/da = dL/db = 0, i.e., 

\J~^gia 7/3 + -^=g 2 ab 2 - Kb = (fJL-1) a, 
3~ 1 

- gi b 7 / 3 + —g 2 ba 2 - na = (jJ,-l)b. 

Eliminating /i from these equations leads to a relation in 
which (a — b) factorizes out, which yields an obvious sym- 
metric solution, a — b. The remaining equation which 
determines the asymmetric solutions is 



9l ab(a 2 l 3 + b 2 ' 3 ) 
x a 4 / 3 + 6 4 / 3 + (ab) 



92_ 

V2 



2/3 _ 



ab — k 



(12) 
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Equations JJjJ and (|12fl constitute a system from which 
a and b should be found, for given gi and g 2 . Although 
this system cannot be solved in an analytical form, the 
critical value, K cr , can be easily found. Indeed, since only 
the symmetric solution exists at k > K cr , the respective 
bifurcation, at which the asymmetric solution merges into 
the symmetric one and disappears, occurs when Eqs. © 
and l|12fl themselves admit the symmetric solution, a = 
b = 7T -1 / 4 . The substitution of this in Eq. (|T2")l yields a 
simple result, 



92 



2gi 



(13) 



15^1/3 ■ 

Note that Eq. 1|13[1 shows that K cr vanishes at the point 



.9 2 = gf = 2 J -^ 1/6 <?i « 0.884 5l . (14) 
V 15 

The meaning of the latter result is that, for g 2 < g 2 , the 
VA predicts the binary fermionic gas to be in the mixed 
state even without the linear-coupling terms. 

Close to the bifurcation which takes place at k = K cr , 
i.e., for < An = K cr — k <C K cr , the system of equations 
© and U3) can be expanded, which yields the corre- 
sponding slightly asymmetric solution, 



a = tt- 1 / 4 ± Aa, b = tt" 1 / 4 t Aa, 
where, in the first approximation, 



(AaY 



Ak 



27V30 



31^32 - 56V2V/6 5l ' 



(15) 



(16) 



It is relevant to note that the condition necessary for the 
symmetry breaking (within the framework of the VA), 
.92 > 9 2 > see Eq. (f 141) . guarantees that Eq. I|16|) yields 
(Aa) > for Ak > 0, as it must be. Farther from the 
bifurcation point, Eqs. © and l|12|) can be solved nu- 
merically, which gives rise to a typical bifurcation picture 
displayed in Fig. where we plot v vs. k for gi = 100 
and gi = 90, 95, and 100. In this case, the fermionic gas 
remains in the mixed state for 172 < <?" ~ 88.4 even for 
k = 0. For 52 > g 2 T , the two species separate, and the 
separation increases with growing 172- Nevertheless, from 
Fig. ^ we nn d that, for g 2 — 90, 95, and 100, parame- 
ter v vanishes for k > 0.68 ,2.67, and 4.67, respectively, 
and the system becomes mixed again. Comparison of 
further analytical predictions with numerical results are 
presented in detail in the next section. 



III. NUMERICAL RESULTS 

We solved coupled MFHD equations Q and J3J) nu- 
merically, using a real-time integration method based on 
the Crank-Nicholson discretization scheme, as elaborated 
in Ref. [3(i|. We discretize the mean- field equations, us- 
ing time step 0.001 and space step 0.05, in domain —25 < 
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FIG. 1: (Color online) A set of bifurcation diagrams — v 
vs. k plots — which show [in terms of the number-of-atoms 
symmetry-breaking parameter, v, see Eq. llOIl l the asymmet- 
ric solution, as found from the variational equations, © and 
(I I 21 . for a fixed gi and varying values of g^. The solution 
with v — is an obvious symmetric one, a — b — 7r~ 1,/4 , see 
text. 



x < 25. We start with the Gaussian ground-state wave 
function of the linear harmonic oscillator as the input at 
t = 0, $(M = 0) = V(x,t = 0) = vr- 1 / 4 cxp(-a; 2 /2), 
while setting g\ = g% = 0. With these initial states, we 
have simulated the time evolution of Eqs. Q and J2J. 
In the course of the evolution, nonlinearity coefficients g± 
and 32 are gradually switched on, at the rate of 0.0001 
per time step of 0.001. Solutions attain the final station- 
ary configuration after the full nonlinearities and linear 
coupling k have been restored. All the wave functions 
displayed below are definitely stable against small per- 
turbations, as stationary solutions to Eqs. Q and 

With k = 0, the norms of the two components are 
conserved separately. In that case, we also used an 
imaginary-time-integration method based on the Crank- 
Nicholson discretization scheme, and the results obtained 
for Eqs. and are in agreement with the integra- 
tion in real time, which verifies the correctness of the 
solutions. However, with k ^ 0, the normalizations of 
the two components are not separately conserved, and 
the imaginary-time-iteration method is not applicable. 

A comparison of the analytical (variational) predic- 
tions from the previous section with numerical results 
can be made in terms of function of g 2 — gf , for 

fixed values of g%. This comparison is displayed in Fig. 
|5Ja), where we plot variational and numerical results for 
gi = 0, 50, 100, and 200. This plot demonstrates that 
the simple VA provides for a qualitatively reasonable pre- 
diction, provided that g 2 — g" and g\ are not too large. 

A simple but noteworthy feature of the VA is the linear 
dependence of g" on gi, see Eq. (TH)) . It can be com- 
pared to its numerically found counterpart, as shown in 
Fig. |2{b). The numerical dependence is nearly linear 
too, with the slope dg 2 j dg\ « 1.3. The comparison with 
the analytical prediction for the slope, which is 0.884 as 
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FIG. 2: (Color online) (a) The dependence of the critical value 
of the linear-coupling strength n CI on the nonlinear-repulsion 
constant <?2, plotted as K cr vs. (gi — g") for gi = 0,50, 100, 
and 200. Numerically obtained curves are labeled by respec- 
tive values of g%, while their counterpart produced by the 
variational approximation variational is labeled "var". Re- 
call g" is the value of gi below which the binary gas remains 
completely mixed in the absence of the linear coupling, (b) 
The numerically found dependence of g^igi), along with the 
respective variational prediction given by Eq. 11411 . The nu- 
merical dependence is (almost exactly) linear, although the 
slope, 1.3, is somewhat different from its variational counter- 
part, 0.884. 



per Eq. I|14fl , shows that the VA yields qualitatively cor- 
rect, but not quite accurate, results. For g\ = 0, the VA 
yields — 0; whereas the corresponding numerical re- 
sult is ~ 7. This means that, even for binary atomic 
systems with no intraspecies interaction (g\ = 0), inter- 
species interaction 172 has to be above the critical value 
g2 ~ 7 for the demixing to take place (which is a conse- 
quence of the fact that the two species can be mixed due 
to the pressure from the trapping potential) . The simple 
VA fails to predict this fact, allowing for the demixing at 
any g% '. 

Besides the conspicuous error in predicting the slope of 
the g" (31) dependence, another shortcoming of the VA 



is seen in the fact that the analytically predicted slope 
of the dependence of K cr vs. <?2 — 9" does not depend on 
gi , while the numerically found slope changes with g\ (in 
this connection, it is relevant to recall that the limit case 
of gi — 0, which is included in Fig. HI (a), corresponds, 
as explained above, not to the MFHD approximation for 
the fermion mixture, but rather to the system of coupled 
GPEs for the binary boson condensate with zero intra- 
species interaction). 

The VA could, in principle, be improved. An essential 
inaccuracy of simple ansatz (jSJ) is the assumption that 
the wave functions monotonically decrease with in 
reality, numerically found wave functions of the ground 
state have several local minima and maxima, see Fig. 
below. Even in the fully mixed state, the profile 
may be markedly different from the Gaussian, as sug- 
gested by the curves pertaining to k = 1.5 in that fig- 
ure. The simplest possibility to improve ansatz JSJ in 
this respect would be to replace it by a generalized one, 
$ = (o + a 2 x 2 ) e-* 2 / 2 , * = (6 + b 2 x 2 ) e'* 2 / 2 , with 
extra constants 02, &2- However, although the Lagrangian 
corresponding to the extended ansatz can be calculated 
in an analytical form, the variational equations following 
from it are so cumbersome that direct numerical solution 
of exact equations (|TJ and J5J), or equivalently, (QJ and 
©, is actually simpler. 

Now we consider the results for the symmetric system 
A/i = 0. Collecting the results for many values of g\ and 
g2, we have concluded that details of the solutions may 
differ considerably, but their qualitative properties are 
universal. An appropriate case to display generic results 
is one with g\ = 142, 32 = 100 (in this case, g" = 137). 
In Fig. |21 we display a set of ground-state wave functions 
found as k increases towards the transition to the full 
mixing at K = k ci (in this ss 1.50). A salient 

(and truly generic, as shown by comparison with a large 
pool of numerical data obtained at other values of g\ and 
(72) feature of the ground state in the immiscible phase is 
that the two species do not collect themselves in different 
spatial domains, separated by a spontaneously formed 
wall (which is the case in the boson mixtures |!4|24(), but, 
instead, they form two spatially modulated anti-locked 
density waves (i.e., ones with a phase shift of 7r), while 
each wave remains an even function of x. In the absence 
of the linear coupling, a qualitatively similar structure of 
the ground state of the binary fermionic gas, with equal 
numbers of atoms in the two components, was reported in 
Ref. |l5j . As k approaches the miscibility point, k ct , the 
modulation in both components becomes more shallow, 
and finally completely disappears for k > K cr , where a 
full mixing between the components takes place. 

As overall characteristics of the immiscible state, we 
use the number-of-atoms symmetry-breaking parameter, 
defined by Eq. J7J, and, in addition, 

^ /- + r* 2 [$ 2 (s)-* 2 (»)] dx 

J^x 2 [<f> 2 {x) + V 2 (x)]dx ' 
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FIG. 3: (Color online). A generic example of a set of ground- 
state wave functions of the binary mixture, Q(x) and ^(x), 
as obtained from the numerical solution of Eqs. Q and 
in the symmetric case (A/i = 0) for gi = 100 and g2 = 142 
and varying k. 



which quantifies the anti-locking between the modulated 
density waves in the two components. For the same typ- 
ical case which is illustrated by the set of ground-state 
profiles in Fig. [21 v and 5 are displayed, as functions of k, 
in Fig. 4. A salient feature of the dependences (which is 
another truly generic feature, observed in all other cases 
for which numerical solutions were obtained) is the oscil- 
latory character of the decay of both v and S, observed 
as k approaches k ci . Naturally, we have v = at k = 
(we always started the computation from a symmetric 
state, with N% = N2, at k = 0, when the two num- 
bers of atoms are conserved separately), while E(k = 0) 
is different from zero, because, as mentioned above, the 
anti-locking between the modulated density waves in the 
two components of the ground state (i.e., demixing in 
this state) takes place also in the system without the lin- 
ear interconversion [l5| . Because the symmetry breaking 
between the two components is spontaneous, the overall 
signs of v and S are, of course, selected arbitrarily. The 
plots in Fig. 4 are generated by choosing S(k = 0) > 0, 
and then drawing the curves by continuity. 

Another feature that deserves special consideration is 



FIG. 4: (Color online) Parameters which characterize the 
spontaneously established inequality in numbers of atoms 
between the linearly coupled immiscible species, v [see Eq. 
@], and anti-locking in the spatial density modulation of the 
species, H [see Eq. 1171 1. are shown as functions of k, found 
from the numerical solution of Eqs. Q and in the sym- 
metric case (A/i = 0). The critical value K cr w 1.50 of k above 
which complete mixing takes place in this case is shown by 
an arrow. 

that, at some values of k (in Fig. 4, these are k w 0.25 
and k ~ 0.75), both v and 5 take nearly zero values 
almost simultaneously, suggesting that the system falls 
into a nearly mixed state at these points. Indeed, the 
consideration of the profiles in a vicinity of one of these 
points, see Fig. shows that the two profiles are almost 
identical at n close to 0.24; nevertheless, the modulation 
does not completely disappear in this case, and the sit- 
uation seems as an accidental mixing. As k moves away 
from this point, both v and 5 become different from zero, 
which implies restoration of the demixing and difference 
in numbers of atoms in the two components, as seen in 
Fig. □ 

It is noteworthy that the deepest anti-locked modula- 
tion in the two components (in other words, a maximum 
of |S(k)|) is attained not at k = 0, but at a finite value 
of the coupling constant. As seen from Figs. Eland 4, in 
the present example this value is k sa 0.5. 



IV. FURTHER RESULTS 

A. The asymmetric system (A/i 7^ 0) 

We now proceed to the role of the chemical-potential 
asymmetry in Eqs. Q and |2"|. If the asymmetry is 
present in Eq. A A* ^ 0, the inequality of the num- 
bers of atoms in the asymmetric components is not spon- 
taneous (i.e., the sign of v is not selected arbitrarily), 
unlike the symmetric system: for A/i > 0, N\ is larger 
than N2 (similar to what was found in the model of the 
asymmetric boson mixture 0). Accordingly, A/i tends to 
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FIG. 5: (Color online). In addition to Fig. [3J an extra set of 
profiles of the wave functions in the symmetric case (A/i = 
0), for gi = 100 and gi = 142, in the ground state around 
K — 0.25, where nearly complete mixing occurs accidentally. 



FIG. 6: (Color online). A generic example of a set of ground- 
state wave functions of the binary mixture, &(x) and ^(x), for 
different k,, as obtained from numerical solution of Eqs. Q 
and in the asymmetric model, with gi = 100, g-z = 142 
and A{i = 0.4. 



suppress the oscillations in dependences v(n) and 5(k), 
which were a salient feature of the symmetric system, 
see Fig. 4. Our results clearly show (Fig. 7) that, in the 
model with a nonzero difference in the chemical poten- 
tials, JVi and N2 are never exactly equal. 

In Fig. El we display a counterpart of Figs. |3J and |3 
i.e., a set of profiles of the wave functions in the ground 
state for g\ = 100 and 92 = 142, in the presence of rela- 
tively strong asymmetry, with A/i = 0.4. Some residual 
oscillations in the shape of the wave functions are still 
seen, but they are much less pronounced than in Figs. |3| 
and[S] 

A counterpart of Fig. 4 for the same example of the 
asymmetric system is displayed in Fig. where we plot v 
and S vs. k for the asymmetric fermion mixture studied 
in Fig. As expected, the figure shows that the oscilla- 
tions in the plots of v{k) and S(/c) are greatly reduced, 
and, on the contrary to the symmetric model (see Fig. 
4), there is no point at which these functions, v{k) and 
S(fc), would vanish; instead, they decay asymptotically 
at large values of n. The nonvanishing of v{k) implies 
that the numbers of atoms in the two components never 
become exactly equal at A/i 7^ 0. 



B. Dynamical relaxation 

The above analysis was based on stationary solutions, 
found from Eqs. and (J5J [the stability of the solu- 
tions was verified within the framework of the full time- 
dependent MFHD equations, and (JS}]. Another as- 
pect of the immiscibility/miscibility transition is relax- 
ation of an initial immiscible state after the coupling 
terms were suddenly switched on. In the context of the 
boson mixture, a similar problem was considered in Ref. 
0. In that work, it was found that the sudden appli- 
cation of the strong linear coupling did not lead to re- 
laxation of the system into a mixed state, but rather 
to persistent nonequilibrium oscillations. Here, we per- 
formed numerical experiments, starting from the phase- 
separated state at k — and suddenly adding the linear 
coupling between the components, with k > K cr (in the 
symmetric system, with A/j, = 0). The result can be ad- 
equately presented in the form of the time dependence 
of the immiscibility parameters, v and H, that were de- 
fined above [see Eqs. Q and i|17fl ]. A typical example of 
this dependence is displayed in Fig. |H] which shows os- 
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FIG. 7: (Color online) The same as in Fig. 4, but for the 
asymmetric system corresponding to Fig. HJ i.e., with A/j, = 
0.4 in Eq. ©. 
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electromagnetic field), which form an immiscible binary 
fermionic gas, may enforce a transition to miscibility. 
The description of the degenerate gas is based on a sys- 
tem of equations for the wave functions of the two compo- 
nents, derived in the MFHD (mean-field- hydrodynamic) 
approximation. The transition to the induced miscibil- 
ity was predicted by means of the variational approx- 
imation, and then verified in numerical computations. 
In drastic contrast with the induced onset of miscibil- 
ity in the bosonic mixture [jj, in the present case the 
immiscible components do not form a domain wall, sep- 
arately filling two domains partitioned by the wall. In- 
stead, they form a pair of anti-locked (7r-phase-shifted) 
spatially modulated density waves, and spontaneously 
break the equality of numbers of atoms in the two com- 
ponents (the latter feature was never observed in the 
bosonic mixture). The relation between the imbalance 
in the spin population, induced by the linear coupling, 
and the spatial density-wave patterns developed by the 
species is similar to the LOFF phase in a bi nary fermion 
gas with unequal values of the Fermi radius (2||- [13] ■ As 
the linear-coupling strength, k, approaches the induced- 
miscibility threshold, K cr , the spatial modulation in the 
density waves vanishes, and the numbers of atoms in 
the components become equal. We have also considered 
the influence of inherent asymmetry between the compo- 
nents (the chemical-potential difference between them), 
and dynamical situations, when k is suddenly switched 
on (in the symmetric system) from zero to a value exceed- 
ing « cr . In the latter case, the system features oscillatory 
relaxation to the mixed state. 



FIG. 8: (Color online) Relaxation of the initial immiscible 
state in the system with gi = 100, <?2 = 140 and re = into 
a mixed state, after re was suddenly switched on from to 4 
(recall that re cr « 1.5 for these values of gi and 52). 



cillatory relaxation to the mixed state characterized by 
v = S = 0. 
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